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Abstract. By considering corrections to the asymptotic scaling functions 
of the photon and electron in quantum electrodynamics with Nf flavours, 
we solve the skeleton Dyson equations at 0(1/Nf) in the large Nf expansion 
at the d-dimensional critical point of the theory and deduce the critical 
exponent /3'(g c ), in arbitrary dimensions, and subsequently present a method 
for computing higher order corrections to f3(g). 
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1 Introduction. 



Recently, various authors have provided an insight into the four loop struc- 
ture of quantum electrodynamics, (QED), by calculating the /3-function to 
this order, [Q]. This renormalization group function, which describes how 
the renormalized coupling constant depends on the renormalization scale, is 
important in relating theoretical predictions of the theory with experiment. 
From a technical point of view, the computation of the functions which arise 
in the renormalization group equation, and also include, for example, the 
anomalous dimension of the fundamental fields, become quite intricate after 
several orders, Indeed the latest calculations could only be completed 
in a reasonable time with the use of computer algebra packages, [0]. To 
push the accuracy of current perturbative calculations beyond the present 
four loop status, an important aid is the provision of independent (ana- 
lytic) analysis. In general, this is not usually possible for all renormalizable 
quantum field theories. For a subclass of theories which possess an internal 
symmetry, however, it is possible to carry out different perturbative analysis, 
which can then be related to the conventional coupling constant perturba- 
tion series. For instance, if a theory possesses an O(N) or U(N) internal 
symmetry, one can let the parameter N become large, whence 1/N becomes 
small, and therefore reorder the perturbation series to analyse all 0(1/N) 
graphs initially, before proceeding to the subsequent order. More specifi- 
cally, this large N expansion corresponds to reordering perturbation theory, 
such that one sums chains of bubble graphs first. As we will be discussing 
QED in detail here, we note that with Nf flavours of electrons this theory 
indeed possesses a large Nf expansion. 

Earlier studies of QED in a large Nf expansion were carried out in H ||] , 
where the authors performed the explicit bubble sum for the photon and then 
used it to renormalize the one loop graphs in minimal subtraction, (MS), 
using dimensional regularization in d = 4 — 2e dimensions. In this way, the 
/^-function as well as the renormalization group function corresponding to 
the electron mass were deduced at 0(1/Nf). Indeed this independent result 
was used in to check the 4-loop term where possible. More recently, 
the electron anomalous dimension was computed in the Landau gauge to 
0(1/N f ) in @ and later to 0(1/N?) in §]. These calculations were carried 
out without having to explicitly sum any bubble graphs and perform the 
usual renormalization. Instead the model was analysed at the cf-dimensional 
critical point where the fields obey asymptotic scaling and are conformal, 
and the theory is finite. Consequently, taking a general power law struc- 
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ture for the fields in the neighbourhood of criticality, where the (critical) 
exponent depends only on Nf and d, and are related to the renormalization 
group function corresponding to the electron anomalous dimension, 7(5), 
it was possible to deduce the critical exponent algebraically by solving the 
truncated skeleton Dyson equation with dressed propagators. Indeed this 
approach was based on the pioneering work of 0, |8) which dealt with the 
O(N) bosonic a model and extensions to other models, H pLQ] . Knowledge of 
these exponents in arbitrary dimensions gives an indirect way of determining 
the renormalization group functions. Since the exponent can be expanded 
in powers of e and the critical coupling g c ~ 3e/Nj, taking a general form 
for say 7(5), its undetermined nth order coefficient can be deduced from the 
coefficient of e n in the expansion of the critical exponent. 

In this paper, we will extend the formalism of || to compute the 
/3-function critical exponent for QED at 0(1/Nf), which is important for 
several reasons. First, it will, of course, provide an independent check on 
the result of H and also ensure that the method we will use is correct. As 
far as we are aware, the proper treatment of the corrections to the asymp- 
totic scaling functions of a U(l) gauge field has not been given previously. 
Having demonstrated its usefulness here, it will open up the possibility of 
carrying out further calculations in other models, such as the bosonic CP{N) 
a model, which possess a U(l) gauge field. Second, and more importantly, 
the formalism we introduce, will provide the foundation for probing the 
/3-function to 0(l/iVj ) and we will outline several technical points in for- 
mulating the algorithm to carry this out. Indeed, we suspect the method of 
Hi, which has proved successful in providing the only 0(1/N 2 ) information 
of critical exponents in other models, will provide a quicker evaluation from 
a technical point of view than §, since, for example, one does not have to 
compute bubble sum corrections with massive electrons. 

The paper is organised as follows. In sect. 2, we introduce the formal- 
ism we will use, discussing in detail the derivation of the asymptotic scaling 
form of the propagators and 2-point functions of the various fields. These are 
used in sect. 3, to represent the Dyson equations at criticality by algebraic 
equations, which can be solved systematically at each order in 1/Nf. After 
pointing out the need to include a higher order two loop graph, which is a 
novel feature for fermionic theories compared with the treatment of bosonic 
theories of |?], ||], we proceed to compute it in sect. 4, where we give an ex- 
tensive discussion about how to calculate certain massless two loop Feynman 
diagrams, using a set of recursion relations. The full calculation is finally 
assembled in sect. 5, where we outline certain technical points concerning 
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the 0(1/ Nj) evaluation of the /3-function. An appendix contains a library 
of the basic results required at certain stages. 



2 Preliminaries. 

We begin by introducing the notation and basic formalism for analysing 
QED with Nf flavours of electrons at its critical point, defined as the non- 
trivial d-dimensional zero of /3(g). First, we recall the (massless) lagrangian 
is 

L = w + A^W - ^gP - (2.1) 

where 1 < i < Nf, is the U(l) gauge field with F^ v = d^A u — d u A^ and 
b is the covariant gauge parameter. We have rescaled the electron charge, e, 
into the quadratic term of the gauge field to ensure there is a unit coupling at 
the 3- vertex, in anticipation of applying the methods of 0, f|. However, we 
will use g = (e/27r) 2 as our coupling constant here, following the conventions 



of [11|. Next, we write down the asymptotic scaling form of the propagators 
of the fields in the critical region, where the theory is massless and has a 
conformal symmetry. In earlier work, we considered only the leading order 
form of these functions which depend on the full dimensions of the fields. 
Here, we go beyond the leading order of || and consider corrections, where, 
as was discussed in other models, [||, the critical exponent involved the slope 
of the /3-function at criticality. Thus, we take the following coordinate space 
forms for the propagators as 

Afj, v (x) r 



B 



2/3 x^Xy 
^ + (2^ -2/3-1) x 2 



+ (x 2 ) x B'(n + 2 ^~ X ) fe^f 

+ 1 ' V» v+ (2/.-2/3 + 2A-1) x 



2 



(2.3) 



The quantities (A,B) and (A',B') are amplitudes which are independent of 
x, whilst a and /3 are the dimensions or critical exponents of the electron, 
■ip 1 , and A^ respectively and 2A = — (3'(g c ), which is a gauge independent 
quantity. 

Several comments concerning (2.3) are in order. First, since we are 
dealing with a gauge field, we have to choose a gauge. As we are carrying 
out a large Nf expansion, this constrains us to use the Landau gauge, because 



4 



from a perturbative point of view renormalization effects are such that one 
remains in this gauge, which is not the case for other covariant gauges. 
Since the large Nf expansion is a reordering of perturbation theory such 
that one sums chains of bubble graphs first, choosing any other covariant 
gauge means one cannot compare results from the large Nf expansion with 
perturbation theory due to the mixing of gauges which would occur, [||, [|] . 
Secondly, with this choice, we note that one can recover the normal Landau 
gauge structure of by transforming (2.3) to momentum space via the 
Fourier transform, 0| ; 

1 _ a(a) f e ikx 



(x 2 ) a 2 2a irV Jk {k 2 y- a 

valid for all a, and its derivatives, where a(a) = T(n — a)/T(a), whence 
(2.3) takes the form 

B 



/i— a 



k^k u 



B' 



(2.5) 



This also illustrates our next point concerning the form of the corrections we 
take for the asymptotic scaling form of the propagator and why the form 
of (2.3) might have appeared contrived. In fact, we have taken the simplest 
correction for A^ u in momentum space having been motivated partly by the 
way one deals with bosonic and fermionic fields. 

From a dimensional analysis of the action, we define, ||, 

a = n + ±rj , (3 = 1-rj-x (2.6) 

where [i = d/2, which introduces the electron anomalous dimension, rj, and 
the anomalous dimension of the 3-vertex of (2.1), \- Both are 0(l/Nf) 
within the large Nf expansion and have been calculated to 0(1/Nf) and 
0(1/Nf), respectively, in the Landau gauge in ||, [7], [DJ, by considering only 
the leading order terms of (2.2) and (2.3), ie 

„ - _ v - _ (2M-l)(2-/i)r(2/i) 

11 Xl ' 4r 2 0u)rGu + i)r(2- M ) [ ' 

which is all we require for this paper, where r\ = J2iZi Wi/Nf. 

A final piece of formalism required for the Dyson equation is the scaling 
forms of the 2-point functions, denoted by tp~ 1 (x) and A~l{x). The formula 
for a fermion field was introduced in as 

~ ^-^ +1 [l-A's(a-l,X)(x^} (2.8) 
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where 



r(a) 



aa(a — /t) 
(fj, — a)a(a)ir 2 ^ 



s(a, A) 



ala 



(a — jjL + X)(a — A) 



(2.9) 



with q(a, A) = a(a — \i + A)a(a — X)/a(a — /x)a(a), and is derived as follows. 
One first maps (2.2) to momentum space and inverts the scaling function 
on this space via G~ 1 G = 1, to give the momentum space form of the 2- 
point function, before converting back to coordinate space, 0. In the case 
of the gauge field one restricts the inversion to the transverse subspace of 



momentum space since this is the only physically relevant piece, [13 
and therefore 



where 



m((3) 



2(2// - (3) x fJi x u 



(2/3 -2^-1) 



(x 2 ) x B'n{f3) Vlu , + 



2(2// - (3 - A) Xfl x v 
(2/3 + 2A - 2fi - 1) x 2 



(2.10) 



m{(3) 
n((3) 



[4(/x-/3) 2 -l]q(/?-/i) 

4vr 2 ^(//-/3) 2 a(/3) 
(/i-/3 + A)(2/i-2/3-l)(2/x 



2/3-2X + 1) 



(/t - (3 - A)(2/x -2(3 + l)(2/t - 2(3 + 2A - t 



q((3,X) (2.11) 



Finally, for completeness we note that the 4-loop (3- function for (2.1) is in 
MS, @, 0, H, 



(d - 4) 5 + 



2N< 



/„2 



f„3 



3- 9 



N f (22N f + 9)_ g4 



64 



616/Vf 2 
243 



+ 



/416C(3) 380\ 
V 9 "27"/ 



144 



I N f + 23 



5 5 + 0( 5 6 ) (2.12) 



Thus, the location of the non-trivial d-dimensional critical point, g c , is 



3e 

w f 



27 e 2 99e 3 77e 4 ^, =, 



4iV/ 



16-/V/ 



16A^ 2 



(2.13) 



3 Consistency equation for A. 

In this section, we discuss the derivation from the Dyson equation of the 
consistency equation whose solution will yield the 0{1/Nf) correction to A 
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and will therefore become the master equation for setting up the algorithm 
for 0(l/Nf) corrections. Moreover, as will be seen, the derivation is not as 
straightforward as in other models. 

First, the skeleton Dyson equations with dressed propagators we initially 
consider are illustrated in fig. 1, which have been truncated to the order in 
large Nf we presently consider. Since we have introduced our formalism of 
the previous section, in coordinate space, one can therefore represent fig. 1 
at the critical point, (2.13), immediately, by replacing the lines of the one 
loop graphs by (2.2) and (2.3), without the need to compute any one loop 
integrals. Thus, 



= r(a-l)[l 
+ 2(2// -l)z 

= A-^x) 



- s(a — 1, X)(x 

(2/z- 2/3-1) 
4zN f 



(3.1) 



[1 + {x 2 ) x A'\ + 



(a?) x B'(J3-n-\ + l) 



(x 2 ) 



2\2a-l 



V/J.U 



(2/z -2/3 + 2A- 1) 
) [l + 2A'(x 2 ) x ] (3.2) 



where z = A 2 B. As we noted earlier, only the transverse part of the Dyson 
equation for in momentum space is physically relevant. So, following 



[13, 14], we isolate this piece of (3.2) by first transforming (3.2) to momentum 
space, noting that we have not cancelled the powers of x 2 , projecting out the 
transverse component before applying the inverse mapping. Equivalently, 
one can make the replacement 



2{a-l)(x 2 ) 



2\a-l 



(3.3) 



valid for all a, in (3.2) before discarding the isotropic tensor and then can- 
celling powers of x 2 . Thus 







m((3) 



(fl-/z) £?'(s 2 )M/?)(/? + A-/x) 
(2/3 - 2/i + 2A - 1) 



(2/3 - 2 M - 1) 



4zK 



f 



(a-1) + A'(x 2 ) x (2a- A - 2) 



(2a - 1 



(2a - A - 1) 



(3.4) 



Next, since both (3.1) and (3.4) contain terms with (x 2 ) x and either A' 
or B' , we treat these portions as independent and decouple each equation 
into two. The terms corresponding to the leading terms of (2.2) and (2.3) 
give two consistency equations whose solution yield r)\. The remaining two 
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equations correspond, naively, to the consistency equation whose solution 
ought to give Ai , by setting the determinant of the 2x2 matrix formed by 
the two remaining equations, to zero. These are 

= [l + „(a-l,AM + ^fc^ (3.5) 
1 V n (2/i - 2/3 + 2A- l)(P-n + l) v ' 

(2a-l)(2a-A-2)A' - (3 - A)(2/i - 2(3 + l)n((3)B> 

(2o-A-l)(a-l) (fx -(3)(2/x- 2/3 -2X + 1) l ' J 

Until this point we have directly followed the formulation used in the 
bosonic and supersymmetric O(N) a models to determine A, [l0|]. How- 
ever, unlike these models there is a fundamental difference in the QED 
formulation. With A = y, - 2 + Y%L\ h/Nj, 

s(a - 1, A) = ^ ~ 2)Nf + O(l) (3.7) 



which is 0(Nf). However, from (3.6) the coefficient of B' is 

y{2fjL — 3)Ai 
"~ 2(/i-l)2( /i -2)(4/ i -7)iV / 



(3. 



which is 0(1/Nf), so that in the final expression one is left with an 0(1) 
quantity involving the unknown Ai. The other elements of the matrix are 
both 0(1). Of course, if one naively follows this path, the correct expression 
for Ai does not emerge, for the simple reason that a contribution has been 
omitted. This arises from the delicate balancing of powers of Nf between 
s(a — 1,A) and the term it multiplies in the determinant to yield an O(l) 
contribution. If we trace the origin of the 0(1/Nf) quantity in (3.6), it arises 
purely from the one loop graph of the equation of fig. 1. However, if one 
considers the higher two loop graph of fig. 2, which would be the next term 
of the Dyson equation for the photon, then it will be of the same order as 
(3.11), and therefore this graph must be included. It corresponds to taking 
the leading term of (2.2) and the (x 2 ) A term of (2.3), in the replacement 
of lines at criticality. For completeness, we remark that this reordering of 
graphs does not occur in the pioneering work of 0, J8|, but does occur in 
other models with fermions, such as the 0{N) Gross Neveu model, though 



the current point was not explicitly discussed in [10]. 

We now amend our equations to include this extra contribution. First, 
we define the value of the relevant two loop correction by 

+ (x 2 ) x B% u (3.9) 
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where Tl^ u = Tlrj^ + Ex^x u /x 2 etc and II, E, ft and E denote the explicit 
values one obtains when the corresponding two loop integral is computed. 
Thus, (3.2) becomes 







AzN t 



'f 



llvK ' (x 2 ) 2 "- 1 
z 2 N f 



2X uXy 



[1 + 2{x 2 ) x A'} 



Ur]^ + E- 



z 2 NfB' 



W\4a+/3-2)i-2-A 



Ur]^ + E 



(3.10) 



Following the procedure we discussed earlier, (3.6) is then replaced by 



8(2a - A - 2) A' 
(2a - A - 1) 

+ zlii 



+ 



2(4a + /3-2/X-A-2 



8n(/3)(a - !)(//- /3- A)(2/i -2/3+1) 
(2a - l)(/i - p)(2n - 2(3 - 2A + 1) 

B' (3.11) 



where it is clear that both contributions to the B' term are the same order, 
0(1/N f ), since ft = S = 0(1) and z = 0(1/N f ). 



4 Computation of 11^. 

This section is devoted to evaluating 11^ by building on our earlier 0(1/N 2 ) 
work, where we developed various integration rules to treat gauge fields 
in the self consistency formalism of ||, as well as to determine some two 
loop graphs similar to 11^. As mentioned, the graph of fig. 2 involves a 
gauge field with an exponent of (3 — A. The shift of the exponent from the 
original has a major consequence. In the absence of such an insertion, 
the graph corresponds to II^, which was discussed in Q, and which is 
in fact infinite due to divergences from vertex subgraphs. However, one 
can compute the divergent piece as well as the finite piece after a suitable 
regulator is introduced. In the present case, the shift to (3 — A means 11^ 
is finite and does not possess the infinities of 11^^, which will be seen from 
the calculation, or by realising that the insertion alters the structure of the 
vertex subgraphs. If, by contrast, one made a shift in an electron line, that 
graph would remain divergent due to the divergence arising from the vertex 
subgraph not affected by such a shift. Whilst II ^ v is finite and therefore does 
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not need to be regularized for subgraph infinites, this means its computation 
cannot be carried out easily. 

In other models, however, it was possible to calculate completely finite 
graphs by the method of uniqueness, which was first introduced in Jl8j and 
developed in various applications later, [19-22]. The basic rule of integration 
is given in fig. 3, where the Greek letter beside each line corresponds to the 
critical exponent for that purely bosonic propagator. When the arbitrary 
exponents, a, (3 and 7, satisfy a + (3 + 7 = 2/i, then one can perform the 
integration over the point of intersection of the three external vertices to 
obtain the product of three propagators with suitably adjusted exponents, 
and represented by a triangle in fig. 3, where v{a\, a.2, 03) = vr^ nf=i a { a i)- 
For QED, it ought therefore, to be possible to calculate the graph of fig. 2 
directly using the same techniques. However, as we noted in ||, there is 
no simple uniqueness rule analogous to fig. 3 for a photon electron vertex. 
Indeed, for such a vertex to be unique, ie integrable, the constraint is that 
the sum of the exponents is 2[i + 2, but from (2.6), 2a + (3 = 2[i + 1 at 
leading order, which is one step away from uniqueness and thus useless for 
our purposes. To calculate ILjj, in ||, in the absence of a suitable rule, meant 
rewriting the integral first by integration by parts rules and then taking the 
trace over the 7-matrices. The former step using the rules of figs. 4 and 
5, reduces the number of 7-matrices or simplifies the propagator, ie the 
first term of the right side of fig. 4. (In figs. 4 and 5, the letter /x or v at the 
vertex of integration denotes the presence of a 7-matrix there, either free or 
contracted with the gauge field. The appearance of a similar letter at the end 
of a line corresponds to a vector z v in the numerator of the integral, where 
Zn is the vector joining this external point with the integration vertex, which 
we represent graphically as a long dashed line.) Having split the integral 
in this way, one takes the trace and then projects out the 77^ and x^x v 
components so that the graph can be written as a sum of purely bosonic 
two loop integrals, where the exponents are related to those of or shifted 
by integers. 

We followed a similar path here to compute II M ^, and if we denote by 
/ the contribution to Yi^ from the first term of the right side of fig. 4 
(omitting the displayed factor), and by the remaining two vertices, we will 
concentrate on ft™. Using the convention trl = 4, together with {7^, 7^} 
= 2if ,v only, then 

IL 1 / = -8{fx-l)[(fi-2)A 3 ^ + B 2 ^ + B 3 ^-2C 3 ^} (4.1) 
II 1 x a xP = -4[S 2 _ /1 + (/i-l)5 3 - M + 2 J D 3 - M + 2E3- M ] (4.2) 
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where we have defined the various linear combinations of purely bosonic 
graphs as 



= 2{a — 1, a, a — 1, a, £) — (a, a, a, a, £ — 1) (4-3) 

B% = (a, a, a, a, £) + 2(a, a — 1, a — 1, a, £) 

— 4(a; — 1, a, a, a, £) + (a, a, a, a, £ — 1) (4-4) 
= (a — 1, a, a, a — 1,£ — 1) — (a, a, a — 1, a — 1, £) (4.5) 

Z?£ = 2(a — 1, a, a — 2, a, £) — (a — 1, a, a — 1, a, £ — 1) 

— (a — 1, a, a — 1, a, £) (4-6) 
= (a — 1, a — 1, a, a, £) — 2(a; — 1, a — 1, a — 1, a, £) (4.7) 

where the general two loop integral (a±, ct2, «3, «4, «s) is defined in fig. 6. 
The reason for these particular linear combinations will become apparent 
later. 

Whilst we pointed out that itself was finite, with respect to vertex 
subgraph divergences, the various integrals making up the definitions in 
(4.3)-(4.7) are in fact divergent when one substitutes the value a = fi, which 
is the leading order value, with respect to Nf, we are interested in. The 
reason for this divergence is that /i is the anti-uniqueness value of a bosonic 
field, and in, say (2.4), we would have terms such as a(fi + rj/2) ~ 
— 2 / (r/T (fx)) . The main point to realise is that whilst individual graphs such 
as (a, a, a, a, 3 — /x), for example, suffer from this infinity, in their sum the 
infinities will cancel and it is therefore important to arrange the calculation 
in such a way that this happens. Indeed this is the reason for the choice of 
combinations in (4.3)- (4. 7). 

We now illustrate this but note that first we need to rewite each integral 
in such a way that the divergence, — a), is displayed explicitly. One 
clue to this is to realise that if each integral could be written in terms 
of integrals with no anti-uniqueness then that graph would itself be finite, 
with, hopefully, a coefficient involving — a). Recently, various recursion 
relations to reduce the exponents of lines in two loop graphs of the form of 



fig. 6 have been given in [22| and used in ||. As they as completely general, 
they can be applied to the present problem and we have illustrated those we 
used in figs. 7-10, where the + or — indicates the addition or subtraction of 
1 to the exponent of the corresponding propagator of the left side. Also, a 
propagator joining the endpoints corresponds to multiplying the integral by 
appropriate powers of x 2 , to ensure the dimensions of each term is consistent. 
To illustrate these points, we derive an expression for Cg by using figs. 
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7 and 8, respectively, for each term of (4.5), to find 

(4a + £ - 2p - 4)(3/x - 4a - f + 3) 



(a -l) 2 
x(a — l,a — l,a — l,a — 1, £) 

e(e + 1 - m) 



(a-l)= 



(a-l,a-l,a- l,a- 1,£ + 1) (4.8) 



where the two terms involving the divergent integral (a — 1, a — 1, a — 1, a, £) 
have cancelled. From (4.8), we see that involves graphs with no anti- 
uniqueness when a = [i and it is therefore finite. For C^-^i, only the second 
term of (4.8) survives and we are left to evaluate (//— 1, /i— 1, /i— 1, /i— 1,4— /i) 
which can be achieved by making the successive transformations /* and \ 
in the notation of |§, to yield (/i — 2, 2, 2, /i — 1, /i — 2). Applying the rule of 
fig. 9 to this leaves the graph (/i — 1, 1, 2, /i — — 1), which has a unique 
vertex. Thus, 

C, - ^ ~ 3)7r2M (4 9) 

In effect the process to evaluate the combinations is two stage. The first 
involves rewriting the constituent integrals to yield finite ones, before eval- 
uating for the specific cases of £ = 2 — /i or 3 — /i, we require. In the 
appendix, we have listed the values of the basic finite integrals analogous to 
{/j,— 1,/i — 1,/i — 1,/i — 1,4 — //}. They were also deduced by the recursion 
relations, which allows one to rewite an unknown integral in terms of others 
in the sequence which have already been evaluated. Further, we have also 
listed the expressions for the remaining combinations, (4.3)-(4.7), though 
we note their derivation was more involved than Cg. 

One graph, (a, a — 1, a — 1, a — 1, £), which occurs in some of the com- 
binations but does not contribute to C^, deserves special mention. Whilst 
it is divergent, we had to use the following circuitous argument to write it 
with explicit l/(/i — a) factors. First, applying the rule of fig. 9, gives 

( a — a — c) 

(a-l,a,a-2,a-l,£+l) + —, r-^ (a-1, a-1, a-1, a-1, £+1) (4.10) 

(a-1) 

Next, we apply the transformation of <— of || to (a — 1, a — 2, a, a — 1, £ + 1), 
use a recursion relation before applying the inverse, <— . Of the three integrals 
which now result, two are finite, (a — 1, a — 2, a — 1, a — 1, £+1) and (a — 1, a — 
3, a — 1, a — 1, £ + 3) with a factor l/(/U — a) in their coefficient, whilst the 
remaining, (a, a — 1, a — 2, a — 1, £+1) still has a hidden divergence. Applying 
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the rule of fig. 10 to this latter integral, the integral we are interested 
in emerges, whilst the remainder are finite. Thus, collecting the various 
contributions and rearranging, we have 

b(a, £)(a — 1, a — 1, a — 1, a, £) 
= c(a, £)(a — 1, a — 1, a — 1, a — 1, £ + 1} 

+ ; r (a — l,a — 3, a — l,a — 1,£ + 2) 

(/i-a) 

(2a — u — 3)<i(a, £) 
+ (2a A - /- m - * - Q { " ~ h " ~ '■ " ' 2 ' " " U + 



where 

6(a,0 



(2a-/i-3)(a-2) 



(2a + £ - - 2)(2a + £ - // - l)(2/x - 2a - £) 



(//-a-Q (2a-//-3)(a-2)(2//-2a-g + l) 
1 (a-1) (a-l)(2a + £-/x-l)(2/x-2a-£) 

v (a - l)(2a + £- /i - 2) v 7 

and we have displayed the l/(/x — a) pieces explicitly. Thus, simply 

(// - lW 2 ^ 

[(/x-a)(a-l,a-l,a-l,a,3-/x)]| a = A( = ^ _ 2)r 2 (,u) ( - 4 ' 13 ' ) 

(it - lW 2 ^ 

[(/x-a)(a-l,a-l,a-l,a,4-//)]| Q = /i = 2 _ ^W^j ( 414 ) 

which completes the basic results required to compute (4.3)-(4.7). Thus 
from the results listed in the appendix, 

fti* = 8^(5^-11^ + 4) 

(M-l)(Ai-2)r2(M) 1 j 

(^-i)(^-2)r 2 (/,) (4 ' lbj 

To complete the calculation, the contribution from W^ v is simplified con- 
siderably by noticing that from the integration by parts rules of figs. 4 and 
5, the factor multiplying IT^, is (fi — a) 2 , so that in applying the recursion 
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relations on its constituent bosonic integrals, very few will survive in the 
limit as a — > /x. Essentially, one need only isolate those terms which lead to 
(/j, — a){a — 1, a — 1, a — 1, a, £). We found, 



~ 7/ 167T 2 " 

"" (M-2)r2(M) 



x 2 



(4.17) 



Thus, overall, 



n = + 

(4/x-7)( M -l)r2( A( ) 1 J 

a _ 9fr 167r^(2/x-3)(/i-2)(/i-3) 

- " w^w^wm ( 9) 

One final point to note is that in the list of finite integrals we used, given in 
the appendix, the function ChT(l, 1) = 3vr 2 ^a(2/i - 2)a(l)[<0'(/i - 1) - ip'{l)] 
appears, where ChT(a, (3) = (a, fj, — 1, fj, — 1, /3, fj, — 1) in the notation of M 
and ip(x) is the logarithmic derivative of the T-function. It turns out that 
this function cancels in the final expression for H^ u . 



5 Discussion. 

Having calculated the graph of fig. 2, we can now evaluate the determinant 
of the matrix formed by (3.5) and (3.11). With the leading values a = fj,, /3 
= 1 and A = /i — 2 for the terms not involving 11^ and noting, ||, 



3)r 2 0u)r/i 

4vr 2 ^(2/x 




(5.1) 



a little algebra gives 



We make several comments. First, it is easy to see that in three dimensions 
the 0(1/Nf) correction vanishes, so that (3'{g c ) = 1 + 0{1/Nj). Second, 
one question which arises in applying the formalism of |7|, || concerns which 
perturbative renormalization scheme the critical exponents relate to. In the 
conventional approach to calculating (3(g), the first few terms of the per- 
turbative series are independent of the renormalization scheme one chooses. 
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Beyond this, at three and higher loops, one must specify the scheme, ie the 
manner in which one removes the divergences to obtain the finte Green's 
functions of the theory. One such scheme is minimal subtraction, where 
only the divergent terms and no finite pieces are absorbed into renormal- 
ization constants. Moreover, this scheme is a mass independent scheme. In 
other models, as well as QED, the exponents one calculates by the method 
we used here, can be related to the corresponding renormalization group 
functions and therefore compared with these functions to the orders they 
are known. So far they are in agreement with MS scheme functions, but 
there remains the question of whether this is true to all orders in pertur- 
bation theory. Intuitively, this is perhaps obvious because the exponents 
are deduced at a point of the theory where the mass vanishes, so that they 
would have to correspond to a mass independent scheme. As the /3-function 
for QED has been explicitly calculated in the large Nf expansion in an MS 
scheme, ||, f|], we can thus compare the exponent we obtained, to that de- 
rived from the renormalization of ||[ and find both are in agreement. 
By the same reasoning, we assume that similar exponents in other models 
calculated using the method of |?], ||] give MS information to all orders in 
perturbation theory. 

The remainder of our discussion is devoted to the algorithm to go beyond 
(5.2), based on the self consistency approach given here. First, the terms 
of (3.5) and (3.11) have to be expanded to the next order in 1/Nf. For the 
one loop graphs and 2-point scaling functions, this is straightforward and, 
moreover, z has been computed to 0(1/Nj) in ||. Next, the correction to 

n^jy has to be calculated, but this is now possible given our discussion of 
sect. 4 where the calculation of the linear combinations, (4.3)-(4.7), were 
discussed in general terms and, therefore, much of the machinery is already 
in place to treat the corrections. The main computing effort lies, though, in 
obtaining the values of the higher order graphs of figs. 11 and 12 and others 
whose evaluation lie beyond the scope of the present paper. These graphs 
fall into two classes. For the two loop graph, we must consider corrections to 
both the electron and photon propagators, (2.2) and (2.3). As we remarked 
earlier, an insertion on the completely internal line yields a finite graph, 
whilst that on a line joining an external vertex will be divergent. However, 
the procedure to regularize and renormalize these divergences was given in 
IJ. For the graphs of fig. 12, which give corrections by the same argument 
that ITyj, was needed for Ai, one need only consider corrections to the photon 
propagator, which means the four loop graph is finite with respect to vertex 
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subgraph divergences. 

We conclude with several remarks. First, we believe that it is now vi- 
able to push calculations in QED to 0{l/Nj) in arbitrary dimensions. We 
have stated earlier why such calculations are required for four dimensional 
perturbation theory, but note that, of course, we will simultaneously de- 
duce useful information about the three dimensional model. Second, as the 
method will be very involved, it might be more appropriate to examine the 
practicalities in a model with similar properties, but with simpler 0(1/Nj) 
structure, such as the O(N) Gross Neveu model, where A2 is unknown too. 
Also, the basic integrals which occur in that model will be relevant for the 
QED calculation. 
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A Various finite integrals. 



In this appendix, we summarize the various building blocks required for Il^ v , 
giving the exact values of various finite integrals first. 



(/x - 1, n - 1, n - 1, n - 1, 4 - fi) 

(n - 2, /x - 1, /x - 1, n - 1, 4 - /x) 
a(4 - m) 



a3(l) 



7T^(/X - l) 2 

vr^r^ - i)r(2 - M ) 



(A.l) 



r(2/x - 3) 



+ (// - 2)'C7iT(l, 1) (A.2) 



(/i - - 1,/x- l,xx - 1,5 -//) = 

(/i - 2,zx - 1,/x - 1,/i - 1,5 - n) = 
{H - 3, n - 1, n - 1, fi - 1, 5 - n) 



tt 2 "(h - If 



a(5 — fx) 



(M - 3) 2 C/iT(l, 1) + 



a(l)a 2 (2) 
(/x - 2, zx - 1, zx - 1, xx - 2, 5 - /x) 



2( /U _3)( /U -2) 2 r 2 (^) 

7T 2 ^-1) 2 

2(/x-2)( A x-3)r 2 (/x) 

7r 2 ^(// 2 -7^ + ll)a(l) 
(/" " 2) 3 a(3 - //) 



a(5 -//)(// -2) 
a 3 (2) 



ChT(l,l) - 



TT^jn - 3)a(l) 
(/i _2)3 a ( 3 - M ) 



(/x - 2, zx - 1, /x - 2, fi - 1, 5 - n) 



a(5 — zx) 



C7iT(l, 1) 



7r 2 ' i (2ii 2 - 12/x + 17)a(2) 



(/x-2) 2 a(3-xx) 
(/i - l,zx - 1,//- l,xx - 1,6 -/x) 

tt^/x- l) 2 ^ 3 -27 / u 2 + 40xz.-24) 
18(/i-2) 2 ^-3)(^-4)r 2 ( M ) 
(/x - 2,/x - 1,/i - 1,6 -/x) 

Tr^/x- l) 2 (3/x 2 - 14/X + 12) 
12( M -2)( M -3) 2 (^-4)r 2 ( / .) 

7r 2 "(/i-4)(/i-l) 2 



(/x - 2, - 1, // - 2, fi - 1, 6 - /x) 
(H - 2, /x - 1, n - 1, fi - 2, 6 - /x) 
(M - 3,// - 1, - 1,6 -/x) 



4(/x 


-2)(/x- 


3) 2 r 2 (/x) 






i) 2 


6(/x 


-3)(/x- 


4)r 2 (xx) 




7T 2 ^ (/X - 


i) 2 


3(/i 


-2)(/x- 


4)r 2 ( M ) 



(A.3) 
(A.4) 

(A.5) 

(A.6) 

(A.7) 

(A.8) 

(A.9) 
(A.10) 
(A.ll) 
(A.12) 
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Next, we give the explicit expressions for each finite linear combination of 
divergent graphs defined in sect. 4, for arbitrary £. 



M = 



2{ijl - a)(2fi -2a- £){2a + £ - p - 1) 



a — 1, a — 1, a, a — 1, £ + 1} 



(a - l) 3 

+ [(2a* - £ - a - l)Gu - a - - (M " «)(« " 1)](3/* - 4a - £ + 2) 
(4a + £ - 2^ - 3) 



(a - 1)< 



■(a — 1, a — 1, a — 1, a — 1,£ + 1) 



(A.13) 



Bp = 



(4a + £ - 2/x - 3)(3/x - 4a - £ + 2) 



(a - 1)' 



(//-!)£(£-// + 1) 
(2a + £-/x-l) 



+ 



+ 2(/i - a 

2(/x - a) 
(a - l) 3 



(a — 1, a — 1, a — 1, a — 1, £ + 1) 

( M -l)£(£- M +l) 



2(// - a)(2a - /i - 1) 



(2a + £- M -l) 



x (a — 1, a — 1, a, a — 1, £ + 1} 



(A.14) 



(4a + i - 2fi - 4)(3/i - 4a - £ + 3) 



(a - l) 2 

— — ^— (a - l,a - l,a - l,a - 1,£ + 1) 

(a — l) 2 



(a — 1, a — 1, a — 1, a — 1, £} 

(A.15) 



Dp 



2(/i-a)(a + g-2) 
(a-l)(a-2) 



(a — l,a — l,a,a — 1, £) 



+ 



+ 



(a-1) 



2 + 



(2a + e-/i-3)(^-a-0 



(a - l)(a - 2) 
x (a — 1, a — 1, a — 1, a — 1, £ + 1} 

2(/i - a)(2/i - 2a - £ + 2)(2a + £ - // - 3) 



(a - l)(a - 2)(2a + £ - A« - 2)(2ji - 2a - £ + 1) 
x (a — 1, a — 1, a, a — 1, £ — 1} 
+ /(a,£)(a-l,a-l,a-l,a-U) (A.16) 



t\ — — ^-J- " '•'« '•'« U • 1 



(a-1) 
2(ji-a) 

Ta^ry 



(a — 1, a — 1, a, a — 1, £) 



(A.17) 
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where we have set 

, . _ (3fi - 4a - g + 3)(// - a - £ + 1)(2// - 2a - g + 2) 
" (a-l) 2 (a-2)(2//-2a-£ + l)(2a + £-//-2) 
x [ ( a _ 1)(2 M - 2a - f + l)(4a + 2£ - 3/x - 4) (2a + £ - 4) 

+ fri-a-£ + l)(2/i-2a-£ + 2) 

x [( M - 0(2a + e - A* - 3) - (4a + 2£ - 3/x - 4)( M - 1)]] (A.18) 

Finally, the finite expressions for the various combinations we needed are, 







2tt 2 ^ 


(A.19) 


A 2 - 




"M = 


(M-i) 2 r 2 ( M ) 


As- 


= 


2vr 2 ^(/x 2 -4^ + 2) 
(/i-l) 2 (/,-2)r 2 (^) 


(A.20) 


B 2 - 


= 


2^(2^ - 3) 
(/.-l)r 2 (/,) 


(A.21) 


B 3 - 


= 


^(2/i 2 + 4) 

(/x- i)(^-2)r 2 Gu) 


(A.22) 


C 3 - 


= 


2tt 2 ^(h - 3) 
(M-2)r 2 (^) 


(A.23) 






2vr 2 ^ 


(A.24) 


D 3 - 


= 


3 " M " r 2 ( / u) 
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Figure Captions. 

Fig. 1. Skeleton Dyson equations for if; 1 and A^. 

Fig. 2. Additional graph contributing to Ai. 

Fig. 3. Uniqueness rule for a bosonic vertex. 

Fig. 4. Integration by parts rule for a gauge vertex. 

Fig. 5. Additional integration by parts rule. 

Fig. 6. General two loop self energy graph. 

Fig. 7. Recursion relation for two loop self energy which increases t\ by 1. 

Fig. 8. Recursion relation for two loop self energy which reduces s\ by 1. 

Fig. 9. Recursion relation for two loop self energy which increases ti by 1. 

Fig. 10. Recursion relation for two loop self energy which reduces t\ by 1. 

Fig. 11. Higher order corrections for electron self energy. 

Fig. 12. Higher order corrections for photon self energy. 
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